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Abstract . The statistics of the energy eigenvalues at the metal-insulator-transition 
of a two-dimensional disordered system with spin-orbit interaction is investigated 
numerically. The critical exponent v is obtained from the finite-size scaling of the 
number Jq which is related to the probability Qn(s) of having n energy levels within 
an interval of width s. In contrast to previous estimates, we find v = 2.32 ± 0.14 
close to the value of the two-dimensional quantum Hall system. 



1. Introduction 

Despite considerable efforts the understanding of the critical behaviour at the disorder 
driven metal- insulator-transition (MIT, Anderson transition) in disordered systems is 
still incomplete. Present day analytical theories are not able to provide quantitative 
results of physical quantities which describe the universal properties in the vicinity of 
the critical point. Therefore, our current knowledge mainly originates from numerical 
investigations using transfer-matrix approaches and Greens-function techniques (see 
e.g. [0). Recently, the energy level statistics has proven to be another powerful 
method to elucidate the peculiar spectral correlations and to calculate the non-trivial 
exponents that govern the non-conventional dynamics H, ^, U, || that are related to 
the multifractal properties of the corresponding eigcnstates. 

In particular, the energy level spacing distribution P(s) is a simple tool to 
distinguish between localized and extended states, and also reflects the respective 
symmetry of the model system under consideration. Here, s = \E — i?'|/A is the 
energy separation of two consecutive eigenvalues E and E' divided by the mean level 
spacing A. On the metallic side of the MIT, random-matrix theory (RMT) serves as 
an adequate description. However, in approaching the transition, novel critical level 
statistics have been found in three-dimensional (3d) |^, ||, |[ ^ and two-dimensional 
(2d) [ij] systems. For small level separations, s, these critical PAs) still 

resemble the RMT- result, P(s) ^ , indicating strong level repulsion, while for large 
spacings a behaviour, P(s) exp(— ks), similar to the uncorrelated Poisson-statistics 
of the localized states is observed. Previously, we have found the size-invariant PJyS) 



for 2d symplectic and QHE systems 1 14 with k k. 4.0, and for 3d orthogonal and 
unitary models where k ~ 1.9 was obtained. 

A disorder driven metal-insulator transition is observed in 3d for each possible 
symmetry class: orthogonal (/? = 1), unitary (/3 — 2) and symplectic (Ji = 4). In 
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non-interacting 2d systems, however, in the presence of both time-reversal and spin- 
rotational symmetry (orthogonal) all states are believed to be localized for any disorder 

> 0. The same holds also in the absence of time- reversal symmetry (unitary), but 
singular energies are present in strong magnetic fields where the localization length 
diverges leading to the quantum Hall effect. 

The only complete MIT in 2d is found for symplectic systems that possess 
time-reversal but no spin-rotational symmetry. Therefore, the symplectic model 
is of particular interest for developing analytical theories and the knowledge of 
characteristic physical quantities is of importance. Yet, there is no consensus about the 
magnitude of the universal critical exponent v that governs e.g. the divergence of the 
localization length, £,{E) ^ [£' — £^o|~'^, at the MIT, whereas the value v « 2.35 for the 
QHE case is commonly accepted |l^ . The proposed values for the symplectic 

system have been found by using the transfer matrix method. They are scattered over 
a wide range: v = 2.05 ± 0.08 at WJV = 5.875 ± 0.010 jl|, |l|l, = 2.75 ± 0.15 
at Wc/V = 5.74 ± 0.03 @, Ip, and = 2.5 ± 0.3 for a different model [||. All 
suggested j/-values are clearly different from the QHE-case which is contrary to the 
observed closeness of the generalised fractal dimensions D{q) in symplectic and QHE- 
systems Since several recent investigations (e.g. |2|, |ll|, |2j]) at the MIT were 
performed taking the numbers {v = 2.75 and critical disorder Wc/V — 5.74) published 
by Fastenrath ^ for granted, an independent check of the critical exponent and 
disorder is necessary. 



Model 



Three different models 26, |1^ have been proposed for numerical studies of the 
localization properties in 2d disordered systems with symplectic symmetry. In the 
present investigation we use the model suggested by Ando in order to be able to 
compare with the above mentioned results for the critical exponent. Also, this model 
seems to be more realistic, because it simulates transfer of electrons between s-orbitals 
via p-orbitals in the presence of spin-orbit interaction. The latter is responsible for 
the symplectic symmetry due to the broken spin-rotational invariance. In second 
quantization, the Hamilton operator on a square lattice with sites m and n and lattice 
constant a is 

H = ^£m4n,a(^m,a+ ^ i^n, a; Ti, a') cl^^^c^y , (1) 

where the disorder potentials are random numbers distributed between —W/2 
and W/2 with probability P{e) — V/W. Periodic boundary conditions are applied 
in both directions. The spin-orbit interaction strength S is defined as the ratio 
S = Vi/{y^ F2^)^/^ where and V2 are matrix elements of the 2x2 complex 
transition matrices V{m, a; n, a'), which depend on the transfer-direction and on spin 
(7. In the following we choose S — 0.5 and the unit of energy V = [V^ + ¥2)^/'^ = 1, 
and only nearest-neighbour transfer is considered. The two-fold degenerate eigenstates 
are calculated by direct diagonalization using a Lanczos algorithm. 



3. Results and Discussion 



The eigenvalue statistics are calculated in a given interval [—0.5 F, 0.5 V] around 
the band centre E/V — 0.0. A careful spectral unfolding procedure was applied to 
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compensate for possible global variations in the density of states which would disturb 
the interesting local correlations. More than 10^ eigenvalues were accumulated for 
each of the parameter pairs {L,W) with linear system size L/a = 15, 20, 25, 30, 35, 
40, 100 and disorder strength ranging from W/V — 5.4 to 6.8 by calculating up to 
3000 realizations of the disorder potentials. As a scaling variable, we consider the 
quantity [ p7| 

/•OO 

Jo{L,W) = / Qo{s)ds, (2) 
Jo 

where Qn{s) for n = is the probability that an energy interval of width s contains 
no energy eigenvalue. It is related to the nearest-neighbour level spacing distribution, 
P{s) = (PQo{s)/ds^ [Q. Compared to previous approaches |^, H, 0] the major 
advantage of our choice is that all calculated neighbouring level spacings are used, 
because no arbitrary cut-off parameter is needed. 

The spectral correlations of a disordered system undergo a crossover from the 
Wigner to the Poisson statistics when increasing the disorder [§[ |[ 0] . In the metallic 
phase, Jo ~ 0.522 is well known for infinite symplectic systems from RMT |2^, and 
in the insulating phase from the Poisson distribution Jq = 1. As expected, the level 
statistic Jo as a function of W changes from the RMT-result {W < Wc) to the Poisson 
limit [W > Wc) continuously for finite L, but discontinuously in the thermodynamic 
limit. It exhibits critical behaviour close to the disorder Wc, which separates the 
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extended and localized regimes. The sign of the finite-size effect is reversed when 
crossing the fixed point, signalizing the delocalization-localization transition. 

Results of Jo{L, W) as a function of disorder are plotted in Fig. 1 for six different 
system sizes. A point of intersection at W/V « 5.98 is clearly observed. Within 
numerical error SW « 0.04, it is consistent with the value reported in but 
certainly larger than that from |2^. We find our Jq{L, VF)-results to fall in between 
the limiting values mentioned above with an size independent J§ = 0.571 ± 0.001 
at the critical point. Similar behaviour is also found for Jn>i{L,W) by calculating 
the probability Qn(s) of the non-zero number of eigenvalues n in the interval s. We 
obtained a scale-invariant sequence of the critical numbers J^, which grow with n, 
Jf w 0.976, J| w 0.996, and converge very fast to unity. In 3d orthogonal systems, 
an analogous critical sequence ||2^ was found earlier. 

Assuming the validity of the one-parameter scaling hypothesis, Jo{L,W) = 
f{L/£^{W)), one can introduce the correlation length £,{W) and re-scale the linear 
size L so that the data of Jo as a function of L/^ for consecutive W overlap with each 
other. In approaching the critical point the precision of ^ decreases due to less overlap. 
Within statistical accuracy all points fall onto one common curve independently of L 
and W. This yields the two-branched scaling function displayed in the inset of Fig. 1, 
where Jo{L, W) is plotted versus L/£^{W). The lower and upper branches indicate the 
extended and localized sides of a complete metal- insulator transition, respectively. 
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Figure 2. The derivative dJo/dW at Wc as a function of the system size L (log- 
log plot). The solid line is a linear fit to the raw data. Dashed and dotted lines 
correspond to the values of the critical exponent v = 2.05 and 2.75 obtained in | ]ie[ ] 
and |20| , respectively. The inset shows the disorder dependence of the correlation 
length ^(1^). 
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In order to calculate the critical exponent v one can expand the finite-size scaling 
function near W = Wc, 

ML, W) « 4 + A{L/aW)y/'' = 4 + A'{W - W,)L'/r (3) 

This gives dJo{L,W)/dW oc L^/r As soon as Wc is detected one can perform the 
two-parameter linear fit on double-log scale, the inverse slope being ly. Fig. 2 shows 
the derivative dJo{L,W)/dW at Wc — 5.98 as a function of the system size L from 
which a critical exponent i' = 2.32 ± 0.14 can be extracted. Recent calculations of 
a symplectic network model yield a similar exponent |2^. As a guide to the eye, 
two additional lines are drawn for comparison, the slopes of which correspond to the 
values of v obtained by Ando |l^ and Fastenrath |l[ , using the transfer-matrix 
method. 

If Wc is not known with sufficient accuracy, a four-parameter fitting procedure to 
Eq. (|^) is required. By applying a ;\;^-criterion similar to that from ^ , we verified 
that the resulting Wc and v do not change considerably. However, their uncertainties 
increase as expected, SW = 0.10 and 6v = 0.16. Thus, our value of v is clearly distinct 
from those obtained previously. Note that the obtained error-bar Sv does not overlap 
with corresponding error-bars of those u from previous studies. The inset of Fig. 2 
displays the correlation length ^ as a function of disorder W. Here, is defined up to 
a constant factor ^o- It diverges at the critical disorder in agreement with the power 
law £,{W) = ^o\W — Wc\^^ ■ It is generally believed that the value of v is universal 
and does not depend on the strength S of the spin-orbit interaction. It remains to be 
checked whether this expectation really holds. 

In conclusion, we have investigated the critical properties of the energy level 
statistics at the metal- insulator transition of a disordered two-dimensional system with 
symplectic symmetry. The scaling function and the critical exponent were calculated 
using the statistics of spacings of neighbouring energy eigenvalues that have been 
obtained numerically. By performing a finite-size scaling analysis we found a critical 
exponent v = 2.32 at a critical disorder Wc/V = 5.98 which is markedly different from 
those values reported previously. Our value of v for the symplectic system is very close 
to that found in 2d QHE-systems, a behaviour that has been observed previously also 
for the generalized fractal dimensions. 
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